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SUMMARY 


A method is described for measuring the dynamic stability derivatives 
of a model airplane in a wind tunnel. The characteristic features of this 
system are that single-degree-of-freedom oscillations were used to obtain 
combinations of rolling, yawing, and pitching motions; that the oscilla- 
tions were excited and controlled by velocity feedback which permitted 
operation under conditions unfavorable for more conventional types of 
oscillatory testing; and that data proqessing was greatly simplified by 
using analog computer elements in the strain-gage circuitry. 

The system described is primarily for measurement of the damping 
derivatives C^p (damping in roll), Cmq + (damping in pitch), 

Unr “ ^np (damping in yaw), and the cross derivatives (rolling 

moment due to yawing) and Cnp (yawing moment due to rolling) . The method 

of testing also permits measurement under oscillatory conditions of the 
static derivatives (rolling moment due to sideslip), Cn^ (yawing 

moment due to sideslip), and Cxn^ (pitching moment due to angle of attack) 

All these derivatives are of particular importance in estimating the short 
period oscillatory motions of a rigid airplane. 

A small number of experimental data are included to illustrate the 
general scope of results obtainable with this system. 

INTRODUCTION 

One important problem in the dynamic motions of airplanes is the 
nature and the stability of the oscillatory modes. In measuring the 
dynamic stability derivatives which apply to these motions there are 
certain advantages in employing oscillation methods in a wind tunnel, and 
the development of such methods has always been attractive to investiga- 
tors. Most of the early measurements of damping in pitch, were made from 
oscillation' tests of a model in a wind tunnel. Damping in roll and 
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damping in yaw have also been measured in this way but, in general, 
experimental difficulties have prevented the wide application of this 
method to the lateral motions. This is particularly true in the case of 
the cross derivatives, yawing moment due to rolling, and rolling moment 
due to yawing, although in one recently developed method (ref. l) the 
yawing moment due to rolling has been successfully measured using a two- 
degree -of -freedom oscillatory technique. 

Most of the studies of the lateral derivatives have been made on the 
basis of steady turning or rolling motions. Theoretical calculations of 
the derivatives are largely based on this assumption. The steady turning 
or rolling flow technique has been used in the systematic studies of the • 
lateral derivatives in the Langley stability tunnel (e.g., refs. 2 and 3 ) • 
Curved or rolling flight is approximated in the test section of this wind 
tunnel by causing the air to follow a curved or spiral path past a fixed 
model. At hi^ speeds, the rolling derivatives have been measvired by 
steadily rotating a sting-mpunted model in a wind tunnel with a dynamome- 
ter and measuring the damping in roll, yawing moment due to rolling, and 
the side force due to rolling. These methods and. other techniques have 
been described and referred to in various NACA publications on the sta- 
bility derivatives for airplane and missile configurations. 


The purpose of this report is to describe an oscillation technique 
-for measuring the lateral and longitudinal dynamic stability derivatives 
in a wind tunnel. It was developed primarily for testing at high sub- 
sonic or supersonic speeds and for this reason three features are believed 
to be of special interest. One of these is the single-degree-of -freedom 
oscillatory system in which various components of pitch, roll, and yaw 
were- obtained by varying the axis of oscillation. Second, the forcing 
system comprised a feedback loop in which velocity feedback was used to 
excite and control the amplitude of the model oscillation. A third feature 
is a system of strain-gage data processing in which electronic analog 
computer elements were used in measuring the amplitude and phase position 
of the oscillatory strain-gage deflections. The advice and assistance of 
the instrument development branch was extremely valuable in developing 
this system of data processing. 

The test apparatus is capable of measuring the moment derivatives 
which arise from angular motions of the airplane. This includes the 
rotary damping derivatives C^p, Cmq + and Cn^. - the cross, 

derivatives Cnp and Ci^ - Ci^; and the displacement derivatives Czp, 

Cnp> and Cmo,. These derivatives are of particular importance in esti- 
mating the short-period oscillatory motions of a rigid airplane. 
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DEFINITIONS AND S'lMBOLS 


Two systems of axes are used in this analysis. The stability system 
of axes with the positive directions of moments and angles referred to 
this system are illustrated in figure 1. The oscillation axes used for 
wind-tunnel measurements are illustrated in figure 2 and defined with 
respect to the stability axes by a set of direction cosines. Primes are 
used with aerodynamic moment and axis designations referred to the oscil- 
lation system of axes . 


The various stability derivatives are defined aS follows: 
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The following symbols are used in the’ report: 



direction cosines between primed and unprimed axes 


moment or product of inertia, depending on subscript, posi- 
tive where negative moment results from positive acceler- 
ation, slug-ft^ 


mechanical spring constant, positive where negative moment 
results from positive deflection, ft-lb/radian 

aerodynamic rolling moment, ft-lb 

aerodynamic pitching moment, ft-lb 
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aerodynamic yawing moment^ ft-lb 

mechanical damping -moment coefficient, positive where nega- 
tive moment results from positive velocity, ft-lh sec/radian 

resistance, ohms 

wing area, sq ft 

torque, ft-lb 

air velocity, ft/sec 

wing span, ft 

mean aerodynamic chord, ft 
voltage, volts 

frequency of sinusoidal oscillation, cps 
transfer function, 

input 

galvanometer current, amp 

strain-gage calibration constant, amp/volt/unit load 

rolling velocity, radians/sec 

pitching velocity, radians/sec 

yawing velocity, radians/sec 

time, sec 

stability system of axes, defined in figure 1 

system of axes used for oscillation tests, defined with 
respect to the stability axes by the direction cosines 

angle of attack, radians 

mean or static angle of attack, deg 

angle of sideslip, radians except where noted 

pitch angle , radians 
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roll angle, radians 
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yaw angle ^ radians 

angle of rotation of model about x* axis, radians except 
where noted 

small angular displacement about y’ or z* axis, radians 

direction angles, defined in figure 2, deg 

phase angles of cr, , and £y* with respect to an arbi- 
trary reference 

air density, slugs/cu ft 

circular frequency of oscillation, 2nf, radians/sec 
wind-off circular frequency of oscillation, radians/sec 

iQ 

dt 

dfn 

dt^ 


( ) maximum value of a sinusoidally oscillating- quantity ( ) 

A( ) incremental value of a quantity ( ) 

Subscripts define the particular axis or motion to which the general 
symbol applies. 


THEORY 

Single-Degree-of -Freedom. Oscillatory System 


The general dynamic motion of a rigid airplane with no moving control 
surfaces requires six differential equations. Three of these define trans 
lation and three define rotation about the center of gravity. If the 
center of gravity of a model airplane is fixed in a wind tunnel, the equa- 
tions involving translation can be eliminated and the motion is defined in 
terms of rotary motions and derivatives by three equations. The system 
can be further restrained so that rotation occurs about one arbitrary 
fixed axis only. In this case the motion is defined by one equation even 
though simultaneous rolling, pitching, and yawing motions may be involved. 

Assume an orthogonal coordinate system, the x*y*z^ system (fig. 2), 
the origin of which is at the center of gravity of the model airplane and 
in which rotation of the model is always about the, x* axis. Equilibrium 
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requires that the summation of the moments about the x'^ y', and z' axes 
be equal to zero. The equation for small angular oscillations about a 
static equilibrium condition can be written in terms of a single variable 

- Ix'x'^* ■ ^x'*^ ” ^x*'^ '^x' ~ ^ (l) 


The sign convention of figure 2 requires that if Ix'x'^ ^x'.> 
are considered positive quantities, their respective moments must be pre- 
fixed by a negative sign since they oppose the motion. The quantity AL' 
is the sum of all aerodynamic moments about the axis of oscillation aris- 
ing from angular deflection, velocity, acceleration, etc., about the 
static equilibrium condition. The aerodynamic moments due to angular 
acceleration and higher-order terms are generally neglected in stability 
calculations, permitting the assumption that 


AL' 


SL’ dL' . 
So "" "■ Sa 


( 2 ) 


= ipv2sb 

and equation (l) could "be written 

Ix'x'’^* + (jx. - t PVSb2 Cz'.) 6 + - I PV^Sb 


a = T^. 


(3) 


(V) 


It is apparent from the left-hand side of equation (4) that Cz'^ is an 
aerodynamic damping coefficient and that a negative value of Ci'. would 

^CT 

result in a positively damped oscillation. A negative value of the coef- 
ficient C z'^ would result in a positive restoring moment about the axis 
of oscillation. This sign convention is thus parallel to that of the 
stability derivatives about the stability axes. 


Equations expressing the equilibrium of moments about the y' and z' 
axes for small oscillations about the x' axis can also be written 


— Ix'y'^ AM* — Kyi€yt — 0 ( 5 ) 

" ^x ' z ' * ” K 2 I € 2 ' ~ 0 . ( 6 ) 

Equations ( 5 ) and ( 6 ) can be written in this simple form only if €yt and 
are sufficiently small compared with a that their effects are negli- 
gible in equations (l), ( 5 ), and ( 6 ) except for the terms Kyteyi and 
. This is accomplished by limiting €y> and £ 2 ' 'to very small 
values but making Kyi and K 21 very large. In other words, the model 
would be relatively easy to deflect about the x* axis but very stiff 
about the y* and z* axes. From a development similar to that of equa- 
tions ( 3 ) and (4) it can be shown that 
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- Ix'y'ff + I PV^Sb C^\d + = Ky-ey. (7) 

” 7x' z ' + ~ PV^Sb ^ = Kz '€ 2* (8) 

The values of the aerodynamic coefficients in equations (1|), (7), 
and (8) will change with the orientation of the oscillation axes in the 
wind tunnel and the attitude of the model with respect to the air stream. 
These changes are related to changes in magnitude and relative contribution 
of the stability derivatives, ordinarily measured about the stability axes 
defined in figure 1. The geometric relation between the oscillation system 
of axes (the x'y'z' system) and the stability axes (the xyz system) is 
completely defined by the direction cosines between the two systems. These 
can be symbolized in the following matrix form 
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( 9 ) 


where, for example, the cosine of the angle between the y and z* axes 
is H. 


The numerical evaluation of these direction cosines is somewhat com- 
plicated since the stability axes do not remain fixed with respect to the 
oscillation axes as the angle of attack is changed, as is apparent from a 
study of figures 1 and 2. It will be shown later in the Appendix that 
certain simplifications are possible in numerical calculations by a less 
direct approach through a set of model axes. Since, however, in the 
present discussion the direction cosines are considered only in symbolic 
form^ it is not necessary to introduce this additional step. 


Small angular motions about the axis of oscillation can be resolved 
into component motions of roll, pitch, and yaw about the stability axes. 
The relative magnitude of each component depends on the direction cosine 


between the x' 
approximations 
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The moments about the stability axes can be expressed in terms of the 
aerodynamic stability derivatives 


AL = i pV^Sb 


It ('=!p<P ♦ 


( 11 ) 


AM = i PV^Sc 


AN = i pV^Sb 


^ (Cm„® + Cjn^a) + Cm^ (Zto.) 




% 


^ (*^np9 + Cn^,iir + Cn^P) + Cn^ .(^) 


( 12 ) 

(13) 


For straight flight, as in the wind tunnel, a = 0 and (3 = -\lr. The aero- 
d;>niainic moments can then he referred* hack . to the oscillation system of 
axes through the direction cosines. 


+ c(am) 

(Ih) 

+ f(an) 

( 15 ) 

+ J(ZM)J) 

( 16 ) 


Thus, the aerodynamic moments indicated in equations (l), (5)^ and 
(6) for oscillation about an arbitrary axis are defined in terms of 
moments about the stability axes by equations (l4), (l5)> and (l6). The 
aerodynamic coefficients which depend on the angular velocity of the model 
can be derived in terms of the stability derivatives as 

c '. = ^ ilL 

PV^Sb A ^ 

2V 


- A^C^p + AC (Cnp + - C^^) + ^ (Cmq + C,j^) + (C^j. - Cn^) 

(17) 


52 


Cm's = ADCZp + CD (Cz^ - Cz^) + EB ^ (C^^ + Cra^) + AFCnp + FC (Cn^. - Cn^ 

(18) 




Cn§ = AGCZp + CG (Cz^ - Cz^) + HB — (Cm^ + Cm^) + AJCnp + CJ (Cnr - Cn^) 

( 19 ) 


Those coefficients which depend on displacement of the model become 
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C » = 2 ^ 

^cr pv^sb 3a 

= - ACCzp + ^ Ch^ 

Cin^y = - BCCjp + EB ^ 

Cn’a = - GCCzp + HB I Cin^ 


- 

(20) 

- FCCnp 

(21) 

- JCCnp 

(22) 


In equation (l7) is the aerodynamic damping coefficient measured about 

the axis of oscillation in the wind tunnel. If the x' axis coincides 
with the X axis the oscillation would be pure roll. In this case 


A^ = 


1 and AC = B^ = = 0 so the measured damping coefficient would be 

, the damping-in-roll coefficient. Similarly, a pure pitching or yawing 


C? 

oscillation would result in the measurement of damping in pitch or damping 
in yaw. 


In general, one stability derivative can be obtained from each separ- 
ate physical measurement. In equations (l?)^ (l8)^ (l9) there are 

eight stability derivatives which depend on angular velocity; however, 
these derivatives form only five independent terms. The derivative 


always appears with C 




^mq 

in the above since, for the pure rotary motions 


considered, q is always equal to d. 
r = - ^ in a test of this type, 


(See ref. U.) Similarly, since 
C appears as one term and 






- as another. The evaluation of these five terms (C 


+ C 




ip' 

C^r - and Cnj. - Cn^) requires five unique measurements. 


Equations (l7), (18), and (19) can be considered in a purely formal 
way as the basis for a system of equations containing the unknown stability 
derivatives. Assuming that five values of or Cj^'^ are available 

from wind-tunnel measurements, along with the appropriate direction cosines 
for the axes about which the measurements were made, a system of five 
equations could be formed. These equations could then be solved simulta- 
neously for the five stability derivatives, providing the equations are 
mathematically determinate . 


The necessary values of or Cnj which lead to the velocity 

derivatives and C^'^, 0^'^, or which lead to the static derivatives 

are obtained from physical measurements of the model oscillation through 
equations (4), (7), and (8). Measurements can be made of the frequency of 
oscillation w, the input torque Tx'. the oscillation amplitude a, and 
the small angular deflections €y' and £2' • There is considerable latitude 
in the choice of axes of oscillation and the particular quantities to be 
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measured within the general confines of mathematical determinateness of 
the stability derivatives. Note, however, that €yi and sire inher- 

ently more difficult to measure than a. The small displacements and high 
stiffness required about the y' and z' a^ces to maintain the validity of 
equations (l^), (j), and ( 8 ) impose a limitation on the accuracy of meas- 
urements about these axes. Friction, backlash, and interaction become of 
increasing importance as the displacement is reduced. Some measurements 
must be made in conjunction with large static pitching moments or aero- 
dynamic disturbances of a random nature and these factors will affect the 
design of the apparatus and the accuracy of the system. 


For the oscillation apparatus discussed herein it was decided to 
avoid measurements of ey< and 6 z' where possible in favor of a. The 
arrangement was such that all the stability derivatives discussed except 
Cn^ and C2^ could be determined from measurements of 10, T^s and 0 in 

a series of four tests in which all' test conditions remained constant 
except for changes in the axis of oscillation. From these measurements 
also, a value could be obtained for the sum Cnp + - Cj*. Therefore, 

one additional measurement was required to resolve these two derivatives. 
This was obtained by measuring €y» during one of the four tests, which 


provided an accurate resolution of C 


Hr 


due to the other aerodynamic derivatives had little effect on €yi 


because in this case the moments 

The 


direction cosines which identify the axes of oscillation used are given 
in the Appendix. These can be used to form the simultaneous equations 
needed to evaluate the stability derivatives. The method used to relate 
the measurements of T^* , o, and Gy» to the derivatives 


etc., is discussed in subsequent sections. 




Feedback Control 


As indicated in the preceding section, measurement of the aerodynamic 
derivatives depends upon an analysis of a single-degree-of -freedom oscil- 
lation defined by equation ( repeated here for convenience. 

Px« - -i PVSb^ cr + (Kx- - i PV^Sb Ci^ 

In the case of a free oscillation Tx* would become zero and the oscil- 
lation would be a damped sinusoid. Use of this method is generally limited 
to test conditions which would not result in oscillatory instability as 
there is no control over the amplitude once the oscillation is initiated. 

For the forced oscillation, Tx» in equation (i]-) is a sinusoidal 
function of time. One case of interest is where the frequency of the 
applied torque corresponds to the undamped natural frequency of the oscil- 
latory system. At this frequency the inertia moments balance the restoring 


^ a = Tx. 


"-X* X 


I a + 
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moments and the final amplitude after the decay of initial transients 
corresponds to a balance between the damping moments and the applied 
torque. The maximum angular velocity of oscillation can be obtained with 
a minimum of input torque at this frequency, as the entire input is used 
to overcome the damping. It is thus a desirable operating point both from 
the standpoint of power requirements and accuracy in measuring the damping. 

One disadvantage pf the forced-oscillation method is that, as with 
the free-oscillation system, testing cannot be conducted where oscillatory 
instability is encountered. At high Mach numbers and high angles of attack 
where minor changes in test conditions may produce large changes in the 
aerodynamic derivatives, a steady-state oscillation is very difficult to 
maintain. In situations such as this, feedback control of the oscillation 
should be considered as it provides a means for automatically stabilizing 
the amplitude and the frequency of the oscillation for any variation of 
damping, either positive or negative, within the capacity of the forcing 
system. 

The system of feedback control used in the present apparatus evolved 
from unsuccessful experiments with the forced-oscillation technique 
described above at high subsonic Mach numbers. After the development of 
the feedback system it was found that Bratt, Raymer, and Miles in England 
had used a similar technique in 19 k 2 but the results of their experiments 
are not generally available. The principle of operation is similar to 
that of the amplitude-stabilized feedback oscillator. 

The oscillatory system was formed by the moment of inertia of the 
model and the stiffness of the restoring springs. Torque was applied to 
this system in the present case through a linkage with an electromagnetic 
shaker. It is convenient to think of the shaker system as a transducer 
which converts an electrical signal input, into a torque. A strain gage 
indicating the angular deflection of the model converted the oscillation 
amplitude into an electrical signal. Feedback was accomplished by using q 
amplified voltage from the strain gage as a source of electrical signal to 
the shaker. ■ Velocity feedback was used in this case and the strain-gage 
signal of oscillation amplitude was differentiated electronically before 
being introduced into the shaker. • ■ 

Thus, for a system with velocity feedback 

Tx' = gCT (23) 

and equation (k) could be written 

Ix'x'cT + ^Px' - ^ PVSb^ Ci^ - g^ CT + ^Kx' - i PV^Sb a .= 0 (24) 

If g and the aerodynamic derivatives are constants, equation (24) is , 
linear. The case of interest is where 

■ Px- - 5 PVSb‘^ c/. - g = 0 


(25) 
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In this case, 


a = ( 26 ) 

The oscillations are sinusoidal and of constant amplitude. The oscillation 
frequency is the undamped natural frequency, given by 


t - i pV^Sb C7' 
2 


The peak amplitude of the oscillation^ a, cannot he defined independently 
of initial conditions if the terms in equation (2^) are constant as 
assumed. Amplitude stabilization would require that the final oscillation 
amplitude be independent of the initial conditions in the same sense that 
a ’’limit cycle” is independent of the starting conditions in a nonlinear 
oscillatory system. The transfer function, g, of the feedback loop can be 
designed to vary with oscillation amplitude in such a way as to produce 
positive feedback at low amplitudes and negative feedback at high ampli- 
tudes with a limit cycle at some intermediate amplitude. This type of 
stabilization, however, would appear to conflict with the requirement that 
g and the other coefficients in equation (24) be constants for sinusoidal 
motion. These conflicting requirements can be satisfied within practical 
limits by allowing g to vary with oscillation amplitude, but at such a 
slow rate that it remains essentially constant through one cycle of opera- 
tion. 



( 27 ) 


A rudimentary circuit of a quasi-linear element which could be 
inserted in the feedback loop to stabilize the amplitude of oscillation 
is shown schematically in figure 3(a) along with a sketch of its transfer 
function, figure 3(t)). The thermister is the nonlinear control element. 

It is characterized by a high negative temperature coefficient of resist- 
ance and as current, either alternating or direct, is passed through it 
the resultant heating causes its resistance to change. The thermal and 
heat-transfer characteristics of the thermister determine the time required 
to reach a new resistance following a change in current. There are many 
variations of the principle illustrated in figure 3 which would produce an 
equivalent result and which can be found in the . literature on amplitude 
stabilization of. electronic oscillators. 

A schematic diagram of the complete feedback control loop is shown in 
figure 4(a) . With this system the input torque, given by gcr, can be made 
equal and opposite to the damping moments acting on the model for any 
value of a by an adjustment of the potentiometer in figure 3(a). For 
amplitudes less than the desired amplitude the damping moment will be less 
than the applied torque and oscillations will build up from rest. For 
amplitudes greater than desired, the damping moments will be greater than 
the applied torque and oscillations will decrease. The only stable oper- 
ating. point is where 
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X' 


i pVSb^ - g = 0 


(28) 


and this can he shown to apply whether the aerodynamic damping is positive 
or negative. 


Analog Computing System 


Using the feedback control system described, the static derivatives 
Cnict^ Cnp, and Cjp can be determined from equation (2?) and an accurate 

measurement of the change in oscillation frequency between the wind-on and 
wind-off test conditions. The equation for can be obtained from 

equation (27) 



-2Kx' 

PV^Sb 



(29) 


Three values of Ci^ are required for different orientations of the axis 

of oscillation. Inserting these values into equation (20) with the appro- 
priate direction cosines provides three equations for the unknowns 



and C 


zp* 


Measurement of the velocity derivatives is more difficult. Early 
attempts to record the output of strain gages with an oscillograph and 
then to measure the amplitude and phase position of each trace proved to 
be an expensive and time-consuming task even with the aid of automatic 
digital computing equipment. The analog computing system discussed 
herein performs the same mathematical processes as the digital computing 
machine, but does so at the time the data are taken and results in a con- 
siderable saving in the time and expense devoted to data processing. 


The measurements required in this case for a determination of the 
velocity derivatives were a, Tjj' , and ey . Each of these time-varying 
quantities can be represented as a Fourier series in wt by the general 
expression 


F(t) = ao + ai cos wt + bi sin wt + . . .an cos nwt + bn sin nwt (30) 


where w is the fundamental frequency of oscillation. Higher-order terms 
are always present to some degree because of buffeting of the model, wind- 
tunnel vibration, etc,; however, only the fundamental component in F(t) 
is of interest. The ’ amplitude and phase position of the fundamental can 
be determined from the Fourier coefficients, defined by 



MCA TN 33^7 


15 


1 

1 F(t) 

cos ujt d(cot) 

(31) 
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« J 

1 F(t) 

sin (ot d(tot) 

(32) 


-Tt 


If strain-gage bridges are located in the oscillation apparatus in 
such a position as to indicate a, Tx*^ and €yt, the output of each gage 
would be proportional to the product of applied voltage and gage deflec- 
tion. Introducing a voltage into each gage of e cos cat results in a 
gage output current of 

i = keF(t) cos cat (33) 

As in equation (30), ca is the fundamental frequency of oscillation so 
that upon expanding^ equation (33) becomes 


i = kefao cos cat + — + — cos 2cat +.bi sin cat cos cat + 
\ 2 2 



A well-damped deflection galvanometer having a time constant much greater 
2tc 

than -jjj- will not respond to currents of fundamental frequency and above. 

Its deflection will be proportional to the average galvanometer current, 
given by 


^av 



(35) 


With equations (34) and (35)^ an expression for ai can be obtained in 
terms of the average galvanometer current 


2iav 


Wherein i^v ^ measured directly at the time of the test and 

k can be obtained from a static calibration of the strain gage. Simi- 
larly^ bi can be measured using a sine wave of voltage in place of a 
cosine wave. The 90^ phase separation between the sine and cosine volt- 
ages was obtained in this case using the input and output, respectively, 
of an electronic integrator. This integrator and other active components 
in the computing circuitry consist essentially of high gain d-c amplifiers 
in which the input and feedback impedances to each amplifier determine its 
specific function. Similar components were used in the feedback loop 
described previously. 
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A schematic diagram of the computing system used is shown in fig- 
ure 4(b). The signal source for the sine and cosine voltages was the 
strain gage, indicating oscillation amplitude, that was used to excite 
the feedback loop. The reversing switch was used to apply the sine and 
cosine voltages alternately to each gage. These voltages were measured 
simultaneously with each reading by means of the rectifier circuits ei 
and eg and the galvanometers. The feedback loop through the attenuator 
was used to suppress any unusually large variations in direct current 
through the integrator, and the capacitors prevented this direct current 
from appeal- ing at the output. 

The in-phase and out-of -phase components, ai and bi, respectively, 
of Tx' , €y' , and a are used to determine the maximum amplitude and rela- 
tive phase position of each. Only the component of T^' and ey< in quad- 
rature with the amplitude is required to calculate the mechanical damping 
and the velocity derivatives. With the notation 

a = a sin (ut + p) 


Tx' = ^x' + v) 

eyi = £yi sin (o)t + |) 


the velocity coefficients for each oscillation condition can be calculated 
from equations (4) and (7) as 


C 


4 

PVSb^ 



T„i sin (v - p) 


0)0 




(37) 


Ky. ey 


PVSb^ 


toa 


sin (| - p) 


(38) 


Four values of and one of were required in this case which, 

with equations (l7) and (l8), yielded the five rotary derivatives Cj , 
*^np) + Cjn^, Ci^ - and Cn^. - On^. 

OPERATION 

Description of Apparatus 

The oscillation mechanism necessary for the dynamic tests was housed 
in a sting- assembly which was matched to the dynamic model and the wind- 
tunnel model support in such a way that it was interchangeable with the 
stings normally used fpr static testing. It was thus possible to measure 
the static force and moment characteristics and the dynamic stability 
derivatives under identical test conditions.. 
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A model airplane mounted on the oscillation equipment in the wind 
tunnel is shown in the photograph, figure 5. An electromagnetic shaker 
was housed in the enlarged portion of the sting downstream of the model 
airplane. Special model construction was required to obtain the necessary 
strength with a minimum of weight since a reduction in weight simplified 
many of the other design problems, particularly those relating to the 
supporting springs. Designed for a wing loading at high Mach nxambers of 
approximately 400 -pounds per square foot, the weight of the model in fig- 
ure 5 is approximately 5 pounds per square foot of wing area. 

A general view of the electronic equipment needed outside the test 
section is shown in the photograph, figure 6. The console on the right 
in the photograph is the power supply for the electromagnetic shaker housed 
in the model supporting sting. The panel rack on the left contains a 
counter for measuring frequency and the various electronic feedback and 
computing elements illustrated in the block diagrams, figure 4. The gal- 
vanometer and read-out system, not shown in figure 6, is the same as that 
normally used for static tests with a strain-gage balance. 

Two oscillation mechanisms were built, one for pure pitching or yaw- 
ing oscillations and one for* combined rolling and pitching or rolling and 
yawing. A close-up view of this latter oscillation mechanism with part 
of the housing removed is shown in the photograph, figure 7* The crossed 
flexure pivots positioned the model and provided the spring restraint for 
the oscillatory system. Two sets of flexure pivots were built which per- 
mitted testing at two frequencies of approximately k and 8 cycles per 
second. A change in the orientation of the axis of oscillation was accom- 
plished by rotating the apparatus about the longitudinal axis of the sting. 
Eight positions were available although only three are unique for the sta- 
bility derivatives considered. 

The strain gage indicating eyf was of the unbonded type since the 
angular deflection about the y’ axis was ±0.0005 radian or less. The 
deflection was held within the above limits by the radial flexures indi- 
cated as the cross-torque restraint in figure 7* This gage was required 
for testing in only one configuration and was mechanically disconnected 
for the other configurations. Because of the extremely stiff cross-torque 
restraint, an interaction correction was found to be necessary. It was 
established from a static calibration that approximately 6 percent of a 
moment about the x’ axis was measured as a moment about the y’ axis 
because of this interaction, and thus a correction to the measured values 

of C was necessary which amounted to approximately 6 percent of the 

m0 

measured values of C7*. . Other interactions were found to be negligible. 

‘■a 

The accuracy. of the data obtained is believed to be sufficient for 
most dynamic stability calculations. Errors directly assignable to the 
computing system are quite small, within 1 percent of the full-scale capac- 
ity and 1° of arc in the phase angle. Other uncertainties were caused by 
variations in the strain-gage calibration constants over a period of time, 
random aerodynamic disturbances especially at the higher angles of attack 
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for pitching or rolling motions, some unavoidable changes in model surface 
condition, mechanical damping corrections, and many other minor effects. 
For example, in a typical case the lateral velocity derivatives Cnp, 

and Cnj. - Cn^ were repeatable to within ±0.005 for a given set 

of conditions, but were subject to a total uncertainty of approximately 
±0.02 as determined from repeated tests over a period of several months. 
Further study and experience with the test apparatus should permit the 
understanding and elimination of this latter type of uncertainty. 


Experimental Data 


Figures 8, and 10 have been prepared to illustrate the general 
scope of data obtainable with the oscillation apparatus described. These 
data were obtained at a low Mach number for the model configuration illus- 
trated in figure 5. Similar data have been obtained at Mach numbers up to 

0 . 95 . 


The effects of oscillation amplitude and reduced frequency at a 
selected angle of attack can be studied from measurements of the type 
shown in figure 8. It is sometimes desirable to measure only the effect 
of frequency or amplitude on certain combinations of derivatives, such as 
Cnp + - C2^ in figure 8, since this can be done with fewer measure- 

ments. These data confirm the assumption of linearity in the small oscil- 
lations of an airplane about an equilibrium position and indicate that the 
effects of frequency on the stability derivatives are negligible for the 
test conditions represented in figure 8. 

Data of the type illustrated in figures 9 and 10 can be used to 
establish the variation of the stability derivatives with angle of attack 
for a mean oscillation amplitude and frequency as this is the form most 
useful in dynamic stability calculations. The data shown in figure 9> 
along with the lift-curve slope, are the aerodynamic parameters of primary 
importance in estimating the short-period dynamic longitudinal stability 
of an aircraft with the control surfaces fixed. The short-period motion 
in this case is essentially a pitching about the center of gravity com- 
bined with vertical translation. The desirability of experimentally sepa- 
rating and and evaluating other derivatives which may affect the 

longitudinal motion depends on the circumstances and on the precision 
required but, in general, the important features of the motion can be 
estimated without these additional aerodynamic parameters. In the stick- < 
free case, a third degree of freedom is introduced by the elevator motion 
about its hinge which may markedly affect the response of the airplane and 
for which the aerodynamic contribution of the free control surface would 
have to be considered (ref. 5) ♦ 

Data of the type illustrated in figure 10 can be used in calculating 
the dynamic lateral stability of an airplane. Analysis of the lateral 
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oscillatory motion with the controls fixed is more complicated than in the 
longitudinal case because of the three degrees of freedom - rolling, yaw- 
ing, and sideslipping. The aerodynamic parameters required, in addition 
to those shown in figures 10(a) and (b), are the side-force coefficients 
due to rolling velocity, yawing velocity, and sideslip (ref. 6). The side 
force due to sideslip can be measured or estimated from steady-flight con- 
siderations alone. Measured values of the side forces due to rolling 
velocity and yawing velocity would be desirable, but in many cases these 
forces are small or can be shown to have negligible effect on the motion 
(ref. 7) • Here again, as in the case of the short-period longitudinal 
motion, free-control surfaces may radically alter the aircraft response 
(ref. 5)- . 

Many of the suggested methods for calculating dynamic lateral stabil- 
ity (e.g., refs. 5, 6^ and 7) do not consider the effects of sideslip 
velocity P because, for typical airplane configurations used in the past, 
these effects have been shown to be small (ref. 3)* This may, however, 
not be the case for current and future airplane types. The effects of 
and Cn^ on the lateral oscillatory motion can be approximated by 

introducing the terms and Cn^ - Cnp into the equations of* 

motion (ref. 6) in place of and This would indicate that, in 

the absence of independent measurements of Cjp and Cnp^ it would be desir- 
able to obtain values of - C^p and C^p - Cnp from oscillation tests 

since this would approximately account for the possible effects of sideslip 
velocity in stability .calculations. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 20 , 1954 
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APPENDIX 


General methods are available for evaluating the direction cosines 
for an arbitrary rotation of one system of axes with respect to another. 
(See, e.g., ref. 8.) In the present case it would be most useful if the 
direction cosines were evaluated in terms of the. angles t] and A illus- 
trated in figure 2. The angle t) represents the mechanical angle by 
which the axis of the crossed flexures is offset from the longitudinal 
axis of the sting and A is determined by keying the oscillation appara- 
tus to the sting in the proper rotational position. The direction cosines 
used in equations (l7)^ (l8), and (20) then become 


A = cos Ojj cos T] - sin sin t] cos A 
B = sin r| sin A 

C =?■- (sin Ojj cos Ti + cos cxjjj sin q cos A) 
D = sin sin A 
E = cos A 


(39) 


F = sin A cos 


In the case of the velocity derivatives^ a considerable simplification 
in the direction cosines can be obtained by referring them to a set of 
model axes which coincide with the stability axes at zero angle of attack. 
The velocity derivatives are then evaluated first about model axes for all 
angles of attack and then transformed to stability axes. Inserting 
c% = 0 in the above expression for the direction cosines results in the 
following values for the tests reported herein, where the double primes 
refer to model axes 



A 

A" 

B" 

C" 

D" 

E" 

pT! 

Type of motion 

450 

0 

1/^ 

0 

-1/v^ 

nc 

t use 

d 

Rolling and yawing 

450 

90° 


iZ-n/T 

0 

0 

0 1 

1 

Rolling and pitching 

450 

180° 

i/s/T 

0 


not used 

Rolling and yawing 

90° 

0 

0 

0 

-1 

not used 

1 1 

Yawing 


Use of these values for the direction cosines resulted in the determina- 
tion of the velocity derivatives about model axes, using equations (17) 
and (18). The transformation from model axes to stability axes was made 
with the following equations where the double-primed coefficients refer 
to model axes. 
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CZp = Cjp cos^o^ + (Cn^ - Cn^) sin^OBi + 


(Cn'p + - Cjp sin Om COS o^n 

CZr - Cz^ = (Cz*; - Czp cos-o^ sin^a^ + 

(Cn'r - - CZp) sin Om cos Om 

Cnp = Cnp cos^om - (Cz'^ - Czp sin^om + 

(Cn'r - Cnp - CZp) sin 0^ cos Om 




(40) 


Cn^ - % = (Cn’; - C^^) cos^o^ + C{^ sin^a^ . 

(Cn^ - Cz'^ + Cz'*) sin c% cos Om 


^mr 






The displacement derivatives ^^3 were not evaluated 

by the above procedure as there was no computational advantage in this 
case. Equations (17) through (22) are developed about stability axes for 
vhich 3 = - -ilf and a = 0. The use of these same relations for the model 
axes system depends on the presence in. equations (17) through (19) of the 
terms due to rolling velocity, and the advantage in using the model axes 
system is that thereby certain terms are eliminated in the equations and 
simple solutions obtained for all angles of attack. On the other hand, 
the use of this relation for the displacement derivatives would require 
the introduction of corresponding terms due to roll deflection about model 
axes. It can be shown that when these terms are introduced, the resulting 
equations are as difficult from the computing standpoint as the direct • 
evaluation of the displacement derivatives about stability axes; therefore, 
in evaluating the displacement derivatives ^^3^ ^^3^ equa- 

tion (20) and the direction cosines for the stability system of axes, 
equation (39)^ were used. 

It is important to note the difference between the model axes system 
used for equation (40) and the system of body axes used in many stability 
calculations. The orientation of the two systems of axes coincides but with 
the body axes system the sideslip angle 3 is defined as the angle between 
the relative wind and the plane of symmetry in the same manner as with the 
stability axes. With the approximations sin o - a, cos a = 1, the side- 
slip angle referred to body axes would become 
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p = - \|f" cos dm + 9" sin c% 

= (- C" cos Otn + A" sin c%) o 

This value for 3 could be inserted in equations (ll), (l3)^ and subse- 
quent equations which would lead to modifications of equations (17) 
through (22) and these new equations would then represent the stability- 
derivatives referred to body axes. Therefore, while there are many simi- 
larities in the two systems, the model axes system used in equation (It-O) 
is not a true system of body axes and should be considered only as a com- 
putational aid. 
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Figure 1.- The stability system of axes is an orthogonal system of axes 
having its origin at the center of gravity and in which the z axis 
is in the plane of symmetry and perpendicular to the relative wind, 
the X axis is in the plane of symmetry and .perpendicular to the z 
axis, and the y axis is perpendicular to the plane of symmetry* 
Arrows indicate the positive directions of motions and moments. 
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Figure 2.- The orientation of the system of axes is obtained by assuming an orthogonal 

system which is originally coincident with the stability system at zero angle of attack and 
with the same positive directions for forces, moments, and motions as being successively 
rotated about the y axis by an angle t] and the x axis by an angle A. 
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(a) Simplified schematic. 


Output 

Input 



(b) Transfer function for several fixed positions of the control potentiometer. 

Figure Simplified characteristics of the amplitude-control circuit. 
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(a) Excitation system. 



Figure 4.- Block diagram of the excitation and computing systems. 
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A-18500.1 

Figure 5.- Model airplane installed on oscillation mechanism in wind- 

tunnel test section. 
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A-19366 


Figure 6.- General view of electronic feedback and computing equipment 

used for the oscillation tests. 
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Figure 7*- Oscillation mechanism with part of the housing removed. 
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Peak oscillation amplitude, a, deg 


Figure 8.- The variation of some of the lateral stability derivatives 
with oscillation amplitude for two values of reduced frequency. 






34 


NACA TN 3347 


OJ 





CQ 

•oa 

O 

I 

u 

a 

o 

*\ 

•CCL 

r«3 

O 


u 

ro 

O 

•V 

P 

u 


ro 

O 


Figure 10.- The variation of the lateral stability derivatives with angle of attack. 
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